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Abstract 



In this work the interaction of electromagnetic field with quasi-periodic media 
has been scrutinized. We have obtained the formula for a distorted medium 



^ \ polarizability tenzor in the X-ray frequency band. Also there have been ob- 

tained the X-rays dynamic diffraction equations for the mediums with arbi- 
trary smooth distortion. For these equations the Takagi-Tuipin equations are 

<—| ■ obtained as a particular case. We have got a very simple formula for the 

coefficient of X-rays reflection on the bent Bragg mirror. We have scrutinized 
the example of calculating the X-rays reflection and focusing by the bent 
Bragg mirror. Also in this work one could find the rezults of calculating the 
I/*} ■ reflected X-rays caustics form. It is shown that there exists a possibility of 

linear dimensions of the X-rays focusing area having typical values of 10~ 7 cm. 
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dia which we define as media obtained by some smooth deformation of perfect periodic 
structures. Mathematically, the term "smooth deformation" means the biunivocal smooth 
mapping of the medium under consideration on the medium with perfect periodic structure. 
^ The main goal of this article is to show that the electrodynamics of quasi-periodic media 

can be reduced to the electrodynamics of the perfect periodic media, and, thus, to make 
possible application to quasi-periodical structures of well-developed classical approaches. The 
correspondence between perfect and quasi-periodical media is possible when the later one is 
considered as the Riemannian manifold with the fundamental metric tensor g a p = 5 a p + 2e a 0, 
where 5 a p is the Kronecker symbol and e a p is the deformation tensor. To clarify this state- 
ment, let us remind, that in the mechanics of continuous medium there are two methods 
describing the deformation processes: the Euler method and the Lagrange method. The co- 
ordinates coincide with the observer's frame of reference in the Euler method. The method 
of Lagrange stands that the co-ordinates are mounted into a medium and are deformed 
together with a medium. It should be emphasized that the co-ordinates of some point of a 



In this article we consider interaction of electromagnetic field with quasi-periodic me- 



medium in the Lagrange co-ordinate system remain unchanged during the medium deforma- 
tion, whereas the co-ordinates of some point of the medium in the Euler co-ordinate system 
are being changed. In particular, the co-ordinates of Bravais lattice points of a medium 
with perfect periodic structure, defined before the deformation by a set of three integers 
n = {ni,n2, ns} will remain constant both during the deformation and after the deforma- 
tion. Hence, the Bravais lattice of a deformed (quasi-periodic) medium will look alike in 
the Lagrange and the Euler co-ordinate systems. However, it is necessary to remember, 
that the area engaged in a quasi-periodic medium in the Lagrange co-ordinate system is a 
Riemannian manifold with the fundamental metric tensor varying from point to point. 

Let us give some examples how the formalism developed in this article can be used 
while describing the interaction of electromagnetic field with media. Generally even perfect 
crystals can be referred rather to the quasi-periodic media, than to the media with perfect 
periodic structure because there are always temperature fields of deformation at least. The 
classic and quantum superlattices electrodynamics of semiconductors also can be referred to 
the electrodynamics of the quasi-periodic media when it is necessary to take into account the 
fields of the deformation, permanently existing in such structures. A big amount of articles 
devoted to the electrodynamics of the carbonic nanotubes is now published. In most cases 
they take into the consideration the perfectly straight nanotubes with the perfect periodic 
structure. Actually there are always some kinds of the nanotube deformation. It is possible 
to reduce the electrodynamics of such nanotubes, and furthermore the electrodynamics 
of ensemble of nanotubes in some matrix completely to the electrodynamics of the quasi- 
periodic media. In real crystals the ensemble of the flaws, dislocations and impurities creates 
a field of deformations, which is usually divided on two parts. The first one is the averaging 
over the ensemble and is the smooth field of deformations. The second one describes the 
fluctuations of the deformations field in comparison with this smooth averaged field. The 
interaction of an electromagnetic field with the crystals, that have the deformation field 
being the averaged smooth field of deformations, can be described by methods suggested in 
this article. 

The last example to be considered is the X-ray diffraction optics of the elastically de- 
formed perfect crystals. Such crystals completely correspond to the definition of quasi- 
periodic media given above. From the point of view of the classical electrodynamics, the 
X-ray diffraction optics is a simple version of the geometrical optics. The only complication 
are coupled waves inside a crystal which appear because of the dynamic diffraction. Despite 
the distorted Bragg mirror is one of the most important parts of almost all X-ray optics 
system, we have not seen papers correctly describing a dynamic X-rays diffraction on elasti- 
cally deformed crystals with mean radius of curvature less or about 10m whereas elastically 
deformed crystals with such media radiuses of curvature are of practical interest. Takagi- 
Taupin equationscJ-B and the series of articles (see bibliography), using this equations as the 
basic ones are valid for elastically deformed crystals with mean radius of curvature more than 
40m, as will be shown in part 4. Furthermore, these papers and also papers, do not show the 
difference between two methods describing the crystal deformation mentioned above (Euler 
and Lagrange methods). The fact is also ignored that the deformed crystal in the Lagrange 
co-ordinate system is the Riemannian manifold with a fundamental metric tensor varying 
from point to point in the crystal. Such disregardings lead to a paradoxical statement, that it 
is impossible to destinguish a perfect crystal and the deformed crystal in the experiments on 



X-ray beam passing through a crystal far away from the Bragg condition. That contradicts 
both daily practice of X-ray beams application and classical electrodynamics because the 
deformed crystal is a nonuniform medium with refraction index and absorption coefficient 
depending on co-ordinates. 

The first section of this article is dedicated to the electrodynamics of media with perfect 
periodic structure and is being auxiliary. As the medium with periodic structure is a medium 
with a spatial dispersion and the relation between the vector of polarization and the electric 
vector is not local, the wave equation is the integro-differential one and is rather difficult 
to work with. In case of media with perfect periodic structure one usually proceeds to the 
(k,u) -represent at ion, where the integro-differential wave equation is replaced by the system 
of the algebraic equations. Then one derives the medium eigenmodes from the dispersion 
equation and gains all key results of the electrodynamics of media with perfect periodic 
structure. Such an algorithm defies the quasi-periodic media generalization, therefore in the 
first section we expound an unknown method, unwieldy and inconvenient on the first sight, 
that allows to substitute the integro-differential wave equation for the equivalent system of 
the differential equations, that describes the electromagnetic field dynamic diffraction by 
the media with perfect periodic structure. This method can be naturally generalized for the 
case of the quasi-periodic media. 

The second section is devoted to the electrodynamics of quasi-periodic media. Proceeding 
to the Lagrange co-ordinate system allows us to apply the formalism of media with perfect 
periodic structure developed in the first section to the quasi-periodic media almost without 
changes in case we consider a quasi-periodic medium in a Lagrange co-ordinate system being 
a Riemannian manifold. 

In the third section we obtain the polarization tensor of a deformed crystal in the X- 
ray frequency band from the first principles by the joint solution of Maxwell equations and 
von Neumann's equation for a density matrix of a crystal. It is necessary to underline 
here, that the dependence of the Debye- Waller factor on the co-ordinates is ignored, the 
Fourier-components of the polarization tensor of a deformed crystal coincide with Fourier- 
components of the polarization tensor of a perfect crystal by the shape only in a Lagrange 
co-ordinate system. These components appear to be very complicated co-ordinate functions 
in the observer system (Euler co-ordinate system). 

In the fourth section we obtain the equations system of the two-wave X-rays diffraction 
on the deformed crystals and after a series of simplifying approximations we show that it 
is possible to receive from these equations a combined Takagi -Taupina-0 equations system 
which appears applicable only for the deformations of the crystals with medial radius of 
curvature more than 40m . 

Finally in the fifth section we use a principle of the geometrical optics locality to offer a 
simple formula for the Bragg reflection coefficient from distorted Bragg mirror and describe 
the X-ray quanta mirror focusing. 

II. PERFECT PERIODIC STRUCTURES 

First we consider interaction of the electromagnetic field with a perfect periodic 
medium. For media with spatial dispersion, the wave equation for electric field vector in 



(^-representation has the form as follows: 

rotrotE(r, u) = ^-E(r, u) + % I d 3 r'x{r, r'; u)E(r', u) + -^j (e) (r, u) 

c J c (2.1) 

where x(v,r';Lu) is the medium polarization tensor, and y e '(v,cu) is the current density of 
external sources. The intrinsic property of perfect periodic media to be translationally 
invariant, makes it possible to expand the polarizability tensor into the Fourier series: 

X(r,r';uj) = x(r + R n , r' + R n ; w) = x(r - r', r'; u) = ^x'(r - r'; u) exp {-i <r (1) , r')} 

i (2.2) 

In the above equations R n = niai + n2a2+ri3a3 is the Bravais lattice vector, n = {m, ri2, n 3 } 
and 1 = {hihih} are integer quantities, t^> = l\h\ + l 2 h 2 + Z 3 b 3 are the reciprocal lattice 
vectors, a l5 a 2 , a 3 and b 1; b 2 , b 3 are the basic vectors of the Bravais and reciprocal lattices, 
corresp ondingly. 

Let an electromagnetic wave be incident on a periodic medium with a given Bravais lat- 
tice. If the wavelength of the incident electromagnetic field is of the order of the elementary 
sell linear extension, solution of Eq. (|2.1| ) is reasonable to seek in the form of the expansion 



E(r, uj) = J2 E i( r > ^)e i(k+r( ' ))r (2.3) 
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where the amplitudes Ei(r, u) vary slowly in the area with the linear dimensions less then 
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the wavelength and (k) = — . By substitution Eqs. ( |2.2|) -( [2~3"D into Eq. (|2.1|) we obtain 

J2 e lT<l)r {L[(k + t®),V\Ei{t, uo) + k 2 ^ ! d 3 r' X l '(r', u)e- l ( k+rW ^E m ,(r - r', u) = 
l i' J (2.4) 

where y e '(r, uj) = and the differentiation operator L[(k + r^), V] is determined by 

L[(k + r (1) ), V] = [V + 1 (k + r (1) )] x [V + % (k + r (1) )] . 

The symbol x stands for the vector product. Expanding the amplitude E 1+1 /(r — r', u) into 
the Taylor series and substituting this series into the integral in expression ( |2.4j ) we obtain 

JdVi 1 '(r', W )e" (k+r(l))r 'E 1+r (r-r' )W ) = x 1 '(k + r(V)E 1+1) (r lW )+ 

t— /(k + rO u) ^^^- --— g-^fk + rW ^ 2 E 1+ ,(r,o;) 

Here and below, the summation over repetitive indices is implicit. Since x a a — 0(l/k 2 ), 
the estimate [Xa,/3\~ l dXa,/3/ dh = 0(1/ k) holds true and, as a result, dx^p/dki = o(x" /3 ). 
Then, in view of that the amplitudes E 1+1 /(r, u) are slowly varying functions, all terms in 
the last expression except the first one can be omitted and, consequently, integro-differential 
equation (|2.4j ) reduces to 



J2 e lTmr {L[(k + r«), V]E,(r, u) + k 2 J2 /( k + A u;)E I+ ,(r, w )} = 0. 
1 i' (2.5) 

Let ro be the radius-vector of the elementary sell center. In that case, any point r inside 
the sell is given by r = ro + rj. Taking this representation into account, let us multiply the 
last equation by exp-f—fr^r} and integrate the product over the elementary sell volume O: 

J2 e *(r(l')-rW)r { f e ^(l')- (1) )^[( k + ^/^ V]E F (r + T], U^f] (2.6) 

v i 

+k 2 J2 / ( k + r(l V) / e^-^E v+f (r + V , cu)d 3 V } = (2.7) 

i n 

In the framework of the slowly varying amplitude approximation taken above, the quantities 
L[(k + t(1 )), V]Ei/(r + T],u>) andEj, j« (r + 7),u) can be treated as constants within the 
elementary sell and, thus, can be factored out from the integrals. In that case, since 
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the equation (|2.5| ) reduces to the system as follows: 

L[(k + r«), V]E,(r, u>) + k 2 £ X ! '(k + r« w)E 1+ ,(r, w) = (2.8) 



Together with Eq. ( |2.3|) , this equation system describes the electromagnetic field dynamic 
diffraction by a perfect periodic medium. 

When the wavelength is much larger then the medium unit cell linear size, electromag- 
netic fields propagating in the medium can be averaged over the elementary sell and, then, 
a hypotatic homogeneous medium with effective constitutive parameters can be introduced 
into consideration instead of the real inhomogeneous one. Polarization of this homogeneous 
medium is related to averaged electromagnetic field by P avg (r,u) = XavgE avg (r, &), where 
Xavg is the polarization tensor which determines the linear electromagnetic response of the 
effective homogeneous medium. In accordance withQ, the effective polarizability tensor 
X av9 (]s.]uj) is expressed in terms of the Fuorier transform of the microscopic polarizability 
tensor (|2.2|) by 



iT g (k;u)=x (k;u) (2.9) 

l + iVx°(k;^) 

The tensor N takes into account interaction between scattering centers and is determined 
by the lattice geometry: 



/g— l>e\T— H, m — rf\ 
d 3 ri'- r'-piv 






lim graddiv / d 7/1 — r (2.10) 

r-,R m to y '| r - R m - rf | 

n 

The function y(r/) is an even function of co-ordinates. In centro-symmetrical lattices it is 
defined by its momenta: 

i |Av^) = 1, i y "dfyjVfa) = 0, ^Jd'rj^rfcpirj) = fi*, • • • . (2.11) 

n n n 

The tensor iV has been evaluated analytically for some types of elementary cellsQ. In other 
cases it can be done numerically. 

III. QUASI-PERIODIC STRUCTURES 

The basic assumption which the procedure of derivation of Eqs. ( |2.8|) -( |2.10|) leaned 



upon, is the relation (|2.3|) accounting for the property of the polarizability tensor of per- 
fect periodic structures to be translationally invariant. Under effect of deformation, this 
invariance, e.g. condition (|2.3|), is disrupted and thus the algorithm of Eqs. ( |2.8|) -([2.9|) 
derivation outlined above becomes invalid. However, an appropriate modification based on 
the approach originated from the continuous medium mechanics can be applied to involve 
deformed periodic media into consideration. In the mechanics of continuous media there 
are two alternative methods of the deformation process description, the Euler method and 
the Lagrange method. In the framework of the first one the co-ordinate system is coupled 
with the viewpoint whereas the second method treats the co-ordinate system frozen into 
the medium and, thus, in the process of deformation the co-ordinate lines are transformed 
along with the medium. As a result, the Euler co-ordinates of a given point x % are changed 
during the deformation whereas the Lagrange co-ordinates of it, £\ remain fixed. Lagrange 
and Euler co-ordinates are related to one another by 

x l = i l + u\i\i\e)\ 2 = 1,2,3, (3.1) 

where w 1 ^ 1 ,^ 2 ,^ 3 ) stand for the deformation field components. In particular, choosing the 
co-ordinate system e, to be related to the elementary translation vectors of the Bravais lattice 
aj, ej = aj/|aj|, we find out the Lagrange co-ordinates of the Bravais lattice sites to be given 
by quantities £ n = \rt\CL\, ^2^2, ^303} both before and after the deformation. Thus, one 
can conclude that in the Lagrange co-ordinates, the Bravais lattice of a deformed periodic 
structure looks like to the Bravais lattice of a perfect periodic structure and consequently, 
the polarizability tensor in the Lagrange co-ordinates proves to be translationally invariant. 
This allows us to state the relations as follows 

x(e,e» = x(e + en,r + en^) = x(e-e , ,e» = ^x 1 (e-e>)ex P {-2<r( i \e / >} 

1 (3.2) 



where r^ 1 ' = {27r£i/ai, 271^2/02, 271^3/03} is the reciprocal lattice covector in the Lagrange 
co-ordinate system; its components in the Euler co-ordinates are defined by Tu(r) = 
(2irlj/aj)dC, :1 /dx l , i = 1,2,3.. It should be emphasized that, mathematically, a deformed 
medium is the Riemann manifold with the fundamental metric tensor as follows 

where £«(£ ,£ 2 ,£ 3 ) is the deformation tensor of the medium: 



1 , n du k du k du l du k 

W + 9jk W + 9lk WW i 



. , . „. 1 . n OU n OU n OU OU , 

eyit 1 ,?,?) = n(9ik^T + 9jk^l + 9ik^^i) (3-3) 



Here g 4 ° = (e,; • e.,) is the metric tensor of elementary cell of the perfect periodic medium. 
Further mathematical analysis accounted for the properties of the Riemann manifolds al- 
lows one to extend the approach being developed to media which need not be obtained 
by deformation of the perfect periodic medium. In order for relation (|3.1| ) to hold true, a 
one-to-one smooth mapping ( |3.2|) of the medium being considered into a perfect periodic 
structure is necessary and sufficient to exist. The media for which such a mapping exists 
we shall refer to as quasi-periodic media. Wave equation for quasi-periodic media in the 
Lagrange co-ordinates takes the form as follows 

+*4 £ M'VWlOlx'K - f; W )<r'^E(<;', w ) + ^£i w «,w) (3.4) 



where g iJ (£) is the tensor inverse to the metric tensor <7ij(£) determined by Eq. 

(9 ij (09 Jk (0 = Sl), <?(£) = det ||^(0||. 

The procedure of derivation of basic equation describing wave diffraction by quasi- 
periodic media is analogous to that applied above under derivation of Eqs. (^.8|)-(P^|). 
As in the case of perfect periodic media, solution of wave equation (|3.4| ) we shell seek in the 
form of expansion analogous to (|2.3|) : 

E(f , u) = J2 E i(^ oj)e^ ki+r ^ i+ikiui ^ (3.5) 

l 

where /c^are the wave vector k components in the Euler co-ordinate system. As before, 
the wavelength A = 2ttc/uj is assumed to be of the order of the elementary cell linear size. It 
should be emphasized that expression ( |3.5| ), as different from Eq. ( |2.3| ), is not the expansion 
in terms of plane waves because the phases $i(k, £) = $i(k, r) = {hi + r i )C + &««*(£) = 

kjX 1 H . t) f l (v) are intricate functions of co-ordinates in both Euler and Lagrange 

dx l 3 

systems. Letting in Eq. ( |3.4| ) y e \^, u) = and substituting Eq.(^.5|) into ([O 
obtain 



£{e^ )+ ^[-e,(0 ® £f (k(0 + r« |)e<(0 + L(k(£) + r« |) + fc 2 ]!^,*,) 

+^ 2 E / Wy/\9(t-?)tf& u^^M-^WV-e^iZ - £', W )} = (3.6) 
i' ^ 

duH£) 
where fej(£) = fe, + % — 7-^ — are the wave vector k components in the Lagrange co-ordinate 

system, fc(£) + r« = (*i(£) + *?>, A*(£) + r 2 (1) ,fc 3 (0 + r 3 (I) ), A = ( * J*., J*_) an d 

differentiation operators L^(k(^) + r^\ — ), L(A;(£)+T^, — ) are determined by the relations 
as follows: 

tfmt) + ^ (I) > f ) = ^'{|j + %(o + ^ (I) ]>^^ + we) + t?]W9 

HHO + r« |) = -^{^ + i[*,(£) + ri l) ]}V99 lJ {-^ + iMO + rf ]} (3.7) 

Integral in Eq. ( |3.6| ) can be estimated in the following way. As has been pointed out 
above, the polarizability tensor x l (k; u>) is proportional to the amplitude of scattering of 
electromagnetic field by elementary cell. By this reason, the quantity x l (k; to) quickly fall 
down with £' and becomes negligible at distances exceeding linear size of the elementary 

(ill 

cell. It means that the Taylor expansion of u\(, — £') = ii*(£) + £, 3 ~Kp + ' ' ' m the vicinity 

of the point £ in the integrand of Eq. ( |3.6| ) can be truncated beyond the first term. Further 
analysis allows one to reduce Eq. (|3.6|) to the system as follows: 

/ ^yi^-OIx 1 '^'; u)e-^^^% +t {e - €, *>) = 



Hence the integro-differential equation (|3.6| ) can be approximated with the good accurency 
by the differential equation 

$>^{[-e,(0 ® A"(k(0 + r« ^)e<(£) + L(k(£) + r« A) + ^] El (£, u) 

+k 2 J2V\9Wx 1 \H0 + T m ;u,)E l+v (Z,u,)} = 
I' 

Just as for the perfect periodic media it can be proved that this equation is equivalent to 
the system 

[-e h (0 ® L\ (k(£) + r« ^)e'(0 + £(k(fl + r« A) + ^] El (£, w) + 

+ ^^i 1 '(kK)+^ ( V)E, + r(^)} = 0. (3.8) 



which describes dynamical diffraction of electromagnetic field by quasi-periodical media. 



mulas (|3.7|) . 

If the wavelength of the incident electromagnetic field is much larger than the linear size 
of the lattice unit cell, the averaging procedure in a quasi-periodic medium in the Lagrange 
co-ordinate system is identical to such procedure for perfect periodic media and leads to the 
equation 



xr&m,<*) = ^ROI%°(k(0;a;) (3.9) 

! + V\9(0\ N (0x °(k'(0;^) 

with the depolarization tensor N(£) determined by 

iV(0=lime,(0®^(k(0 + r(l),^) 

r , e -<*ie-e'+u(o-u(€ , 4i;)i 

x y ^wes ^ ^_ f + u(0 _ u(e + g _ - r ^)e'(o « (3.10) 



lim e,(0 ® Lf(k(0 + r(l), A) f <Pr,y/\g^ + r,)\ ^Wt n T e ^ 

where the integration is over unit cell of the perfect periodic medium and the function f{rj) 



is determined by the equations (|2.11| ). The equations (|3.9| -??) together with the equation 



W0» ^J^|r^e'(0 - -L^' JriEK,.) = 

^EK,^) + ^r" 9 «,k(0lw)EK,^)+^j«»>K,^) (3.11) 

c c c 

are determining the electrodynamics of the quasi-periodic medium in this frequency band. 



IV. THE POLARIZABILITY TENSOR OF ELASTICALLY DEFORMED 
CRYSTALS IN THE X-RAY RANGE 

In previous section we have obtained the basic equations of electrodynamics of quasi- 
periodic media. Now, starting with the Maxwell equations and von Neumann equation for 
the density matrix, we derive an explicit expression for the polarizability tensor in deformed 
crystals, which constitute a significant class of quasi-periodic media. Since the crystal unit 
cell is about a few Angstroms in size, dynamical diffraction by crystals manifests itself in the 
X-ray frequency range. Below we restrict consideration to this frequency range. As has been 
shown above, a deformed crystal in the Lagrange co-ordinate system is similar to a perfect 
crystal in the Euler system; consequently, the Lagrange co-ordinate system is preferable for 
evaluation of the polarizability tensor of deformed crystals. 



The current density induced in the point £ = (£\£ 2 ,£ 3 ) by X-ray quanta in elec- 
tron subsystem of the crystal can be determined by the equation j^(£,t) = ( j '(£,£)) = 
Sp(x(£, t)p(t)), where (3 = 1,2,3, p(t) is the density matrix of the crystal, J (£,,t) = 
2.3 (£ — £sj t) is the /3-th component of the current density operator in electron subsystem 

s 

(summation is carried out over all crystal electrons). In the nonrelativistic approximation 
for the Dirac equation, the current density operator of the electron j (£ — £ s , t) is given bya: 

fit - &, *) = -— n(e - 6)^(e, t) + fa - o, 

771C 

/(£ - 6) = e[<A + ^V M ^ 7 ][<K£ - ^ + |^(£ - £ s )] (4.1) 

e 2 

Here the value ??.(£— £ s )/r(£, t) is the potential (Rayleigh) part of the current operator, 

mc 

n(£ — £ s ) = 5(£ — £ s ) is the density operator of the s-th electron, A^(£,t) is the /3- 
th component of the electromagnetic field vector potential, j@ (£ — £ s ) is the sum of current 
operators of resonance electric and magnetic transitions, e^ is the asymmetric tensor (£123 = 
I); Pv = —iK^n is the momentum operator in the Lagrange co-ordinates system (V M is the 
covariant derivative), and a 7 , 7 = 1,2,3 are the Pauli matrixes. 

In the interaction representation, the von Neumann equation for the density matrix of 
the crystal interacting with X-ray quanta has the form as follows: 



am 

at 

H'(t) = e^^ T (t)e-t^. (4.2) 



ih^ = [H' ei (t),p(t)] 



Here H stands for the unperturbed Hamiltonian of the crystal while H ei (t) describes its 
interaction with X-ray quanta, 



H e7 (t) = -iyVev^fc&tM'&t)- 



Accordingly to the perturbation theory the von Neumann equation solution can be repre- 

00 

sented by the series pit) = } J P jt), where p^ k \t) is proportional to the fc-th power of the 



fc=0 



electromagnetic field: p (fe) (t) ~ (H' ei (t)) k ~ (E) k . Substituting this series into Eq. (fO 
and equating the terms of the same power in H' (t), we obtain 



Note that H' n {t) = and — ^— = in the absence of X-ray quanta, i.e., p(t) = p , where 



: dpW(t) 
~dt 

dp{t) 

dt 
Po is the unperturbed crystal density matrix. Further p is assumed to be given by the 

F-H Q 

equilibrium (Gibbs) density matrix, p = e kT , where F is the free energy and T is the 
crystal temperature. The average value of the current density operator in the interaction 
representation can be obtained from the formula 

10 



/(£,i) = Sp(ei™f(U)e-^ ot p(t)) = Sp(f H (U)Po), 
where f H {tl,t) is the current operator in the Heisenberg representation: 

OO t t\ tn-1 

+ E (^ / ^ j dt * ■ ■ ■ J ttnl- ■ ■ le** ot ?& *)e-* Aot , «*i)L • • • , ^ 7 (^)]. 



n=0 



Here the time moments are put in the order: t n < t ra _i < • • • < t\ < £. 

In the linear approximation with respect to electromagnetic field, this equation reduces 
to: 



e 2 



/(£, t) = J2 Sp(e f " ot n(Z - tie-WMAPfo t) 

mc z — ' 

s 

t 

S 8 -OO 

SpUe - ZM^-^f{Z - We-i^^lM"^ *i), 

where the summation is carried out over all crystal electrons. Let us present co-ordinates 
of the s-th electron in the following way: £ a = £ n + 6 + £fc + ^( n ) where £ n stands for the 
co-ordinates of the unit cell which the s-th electron belongs to, £/ are the co-ordinates of the 
/-th atom inside this unit cell. In such situation, the identity \. = /, /, /, holds true 

s n I k 

and, thus, the density of current induced by the X-ray quanta in the crystal is represented 
by 

/(£, t) = -— J2 E Sp(ei A °%(t - & - 6 - 77(11, Q) e -**>«po)^& *) + f E E 

n Z n i 

oo 

| dt' j d 3 C'V9¥){Sp[e^ 0t 'j?(C ~ U ~ 6 " »7(n, OK** * ji,^ - & - 6 - »tfn, 0)Po] 
o 

-Sp[M? - e- " 6 " ^(n, 0)e* V jf « - £n - 6 - ^(n, OK^'po]}^', t - f )■ 

Here hi(£ - £ n - & - 77(11, /)) = ^ n(g - g n - 6 - & - v( n J)) is tne density operator 

fe 
of the electrons in /-th atom and jf(£ — £ n — ^ — 77(11,/)) = /_]/(£ — £ n — 6 — £fc — 

fe 
77(11,/)) is the current density operator in /-th atom. In this equation we have neglected 

the components Sp < [jf (n, t'), j/' iM (n', 0)]po | , n ^ n' or I ^ f describing the nonradiative 

migration of excitation from /-th atom to /'-th one. In another words, we have neglected 
derealization of exciton in crystal. Besides, we leave out the peculiarities of X-ray quanta 
interaction with the band electrons, assuming all the electrons being localized on the atoms. 
These approximations are proved to be correct for situation when the incident photon energy 

11 



exceeds significantly the energy of any resonance transition in the crystal. Indeed, this 
condition allows us to restrict ourselves to the impulse approximation^ where there is 
no difference between energy spectra of band electrons and electrons localized nearby the 
nuclei: the electrons density periodicity in the crystal is of the only importance. That is 
why, describing the interaction of X-ray quanta with the crystal, we use the terminology 
appropriate rather to the isolated periodically arranged atoms devoid of zone structure then 
to the crystal. It should be noted that the above assumptions are not applicable to the 
X-ray spectroscopy where both electron derealization and peculiarities of the X-ray quanta 
interaction with band electrons should be taken into consideration. However, this problem 
is beyond the scope of this article and will be considered elsewhere. 

The remarks made above allow us to present the crystal Hamiltonian H by H = 
y y Hi + H p h, where H P h is the Hamiltonian of the phonon subsystem and Hi is the 

n I 

Hamiltonian of the /-th atom in the unit cell. 

For the atomic ground state in the crystal it is possible to neglect the influence on 
this ground state of elementary crystal excitations (including influence of interaction with 
phonons ) . In such a case, the crystal density matrix p can be presented by the product p — 

Pq pQ of electron subsystem density matrix /3q = e n ; and phonon subsystem 

density matrix /3q = e^' ~ ph \ where (3 = 1/kT, and F^ e \ F^ p ' are the free electron and 
phonon energies. As different from the ground state, the contribution of crystal elementary 
excitations can not be neglected for excited atoms because they lead to the appearance of 
the new decay channels of the exited state and, consequently, increase the T-width of the 
exited state. Therefore, we will assume the Hamiltonians of ground and exited states to be 
different. 

Now, let us discuss how the phonon spectrum is modified under effect of the crystal 
deformation. In the harmonic approximation, the phonon subsystem Hamiltonian of the 
deformed crystal in the Lagrange co-ordinate system takes the form as follows: 

H Ph = E 5> Q/3 (n, P " (n,/ ]f (n ' + E E $ ^ n > n > *> 'V(n, W, I')} 

n I l n' I' (4.3) 

where p a (n, I) is the a-th component of the momentum covector of the I -th atom in the n- 
th cell presented in the Lagrange basis. Quantum-mechanically, the momentum is given by 
p a (n,l) = —ihV a with V a as the covariant derivative. Further we assume the deformation 
tensor to be constant over the unit cell; then, g al3 (n,l) = <5 a/3 + e a/3 (n, I) w g a/3 (n) . This 
approximation allows us to reduce the Hamilton canonical equations to 

tu 2 B«(n, I) = <ft>{n) E E D "^ ^ ; > Z ') S >'' ''), 

n' V 

where B a (n, I) = y/Mrfin, l)e 1 ^ and D aP (n, n'; I, I') = ^^'".V'^ = Atf(n - n'; I, I') 

V MiMy 

is the dynamical matrix of the crystal given in the Lagrange co-ordinates. The above 
equation for B a (n, I) can be rewritten as 
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U> 



''B a (n,l) = (^ + ^)^^D Q/3 (n-nV,/ / )5 /3 (n , ,0, 



n' l' 



where the deformation tensor e al3 is considered as a parameter. Then the translation invari- 
ance of Hamiltonian Q allows us to write down B^(n', f ; e) = e ~*e (*-*') _B«( n) I'; e); this 
leads us to the equation 



<jj 



! B a (n, 1; i) = (6 a i + e a ^) £ D 7/3 (q; I, l')B^n, I 1 ; i), 



where _D 7/ g(q; /, /') = } ^-D 7/ 3(n; /, l')e l *"w j s the crystal lattice dynamical matrix in the 

n 

Fourier representation. Hence, we immediately obtain the dispersion equation 



det 



[0 



!rt M , _ $«i J2 D^(q; I, I') - e^ ]T D 7/3 (q; /, t) 



0. 



(4.4) 



As the Bravais lattice of the deformed crystal in the Lagrange basis coincides with the Bra 
vais lattice of the perfect crystal in the Euler co-ordinates system, the matrix ||Z) 7/ g(q; /, Z')| 
coincides with the dynamical matrix of the perfect crystal. Consequently, dispersion equa- 
tion ( |4.4| ) does not contain undefined quantities. Since the matrix (5° 7 + £ a7 )-D 7j g(q; /, I ) 

is Hermitian, the equation roots u>j(q, i), are real quantities; here j — 1, • • • , 3r with r as 
the number of atoms in the crystal cell. 

This matrix is Hermitian, its eigenvectors corresponding to these roots satisfy the con 
ditions of orthonormalization and fullness 

Y,e* a (q,]]l;i)e p (q,]-J';e) = 6%5 hl , 
j 

^2e* a (q,j;l]i)e a {q,j']l;i) = 5 jd > 
i 

i 
and allow us to implement the normal modes of the crystal e a (q,j;n,l;i) = 

(q, j; I; e)e lki ^^ n > . For these modes we have 

Yl e * a (^' 3\ n 5 h e)e/s(q, j; n', l'\ i) = 8p5i,i>5 n , n > 
id 

Yl ^fa^'i n > h s)e a {q\j'\ n, I; e) = <^A(q - q') 



N 



n.l 



where A(q) = / j ^q,r(i) an d r(l) = ( -. -. ) is the vector of the crystal reciprocal 



5 5 

Cl\ 0,2 O3 



lattice. Next, using the standard method we obtain the Hamiltonian of the crystal phonon 
subsystem in the harmonic approximation 
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H ph (e) = J2^M^)[bUi)KA^) + ~] (4-5) 

id 

where b • (e) and 6 qj - (e) are the phonon creation and annihilation operators. It is easy 
to verify that the proposed algorithm of accounting the influence of deformation on the 
crystal phonon spectrum also suits the case of anharmonic oscillations of the crystal lattice. 
At the end we will make two remarks. The first: since the deformation tensor e a ^ = 
e al3 {£) is the essence of the smooth co-ordinates function, the equations (|4.4j)-(|4~5|) are local. 
The discussed above method of accounting of deformation influence on the crystal phonon 
spectrum is the approximation, that can be considered precise enough only when the length 
of the phonon free path in the crystal is less then the linear dimensions of the area where 
the deformation tensor variation can be observed. The second remark concerns the fact that 
the introduced parametric method of accounting the deformation influence on the crystal 
phonon spectrum is not something especially new in the solid-state physics: this method is 
being the standard one in accounting the deformation influence on the energy spectrum of 
band electrons (see for example^). 

After the discussion made above let us continue the calculations of the current in the 
deformed crystals that was induced by X-ray quanta. First let us take up the Rayleigh 
component of the current 

Jm(^ t) = -— J2 E S P( e ^ 0t ^ " & - 6 " vfa J))e-**°*A>)^& t), 

rnr < ' < ' 



mc 



luJ 



that makes the main contribution. As e n Hot p e hHot = po and E(£, ui) = — A(£, u>) 

c 
in the Coulomb calibration, using the Fourier series hi(£ — £ n — £/ — 77(11, Z)) 



/AWk^.^.^^M, we obt ai„ the formula for Rayleigh component of 

he current 

&&<*) = — /o/ 3 ke^X e ~* e(n) £ e ~ J ^ 



The value Sp(e- ikivi{n ' l) p^ h) ) = e - Wl( - k '° is the essence of Debye- Waller factor where W t (k, £) 
is to be determined by the equation 

wi(k, = -A_ y \^(y,hi)\ 2 coth MM (4 6) 

qj J 

Accordingly the Debye- Waller factor in the deformed crystals is the function of the defor- 
mation tensor and consequently is the function of co-ordinates determined by the equations 
(|4.4|)-( ^3|) . As Sp(hi(k)pQ ) = //(k) is the form factor of the l-th atom unit cell, the stan- 
dard formula %, e_ * — N ce u N^ <5(k — r*- '), where N ce u is the number of the unit cells 

n 1 

in the volume unit, allows us to obtain the final formula for the Rayleigh component of the 
current 
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moo *-^ *-^ 
1 / 

or the formula for the Rayleigh component of the deformed crystal polarization tensor in 
the X-ray frequency range 

'"■"" xi (4.7) 



mw z 



In the frequency range where the X-ray quanta energy is greater then the binding energy 
of any electron in the crystal, the current component describing the absorption of quanta by 
the atoms and the consequent decay of the exited state which atom passed into by means 
of quantum absorption is usually much smaller then the Rayleigh component (in the X-ray 
optics it is called the dispersion correction). Nevertheless the accounting of this component 
is very important as it determines the absorption of X-ray quanta in this crystal. Using the 

Fourier decomposition jf(£ - £ n - £, - 77(11,/)) = [ d 3 kfl (k)e^V^ i(n) e-^e-^ i(n,l) , 

this current component can be written down in the following way 



00 
&p& ") = ^ / d ^' V^O E E / dt> J d " k J Me^e^'e-tto+WWe 



-Hki+k'je, 



n I Q 



In the X-ray frequency range under the consideration T (the width of the exited atom 

state in the crystal) is about 10 16 — , that is why in the X-ray diffraction optics the as- 
sec 

sumption Hi + H p h(^) w Hi is being considered. This assumption means that we neglect 
the influence of Raman scattering of X-ray quanta on the phonons. If we take this as- 
sumption and use the standard formulas } j e~ l ^ ki+k ^ ( n > = N ce u Nj(k + k' — r^ 1 ') and 

n 1 

Sp{e- i{ki+k > i{n ' l) pi h \i)) = e -^ (k+k '' §) , we obtain for jjL,(£,w) the following equation 



JdispV* ' 



U> 



fojj 

Ja(l) 



^ e -^(r«,0 Y^e-^M J^./^e-ir^ J^^^-^)_^_ 



J2 £>(A(a))[<a(0; J«(0,M,|jf(k) / ^ e -i^-^-^'i^(-k-T«)|a(0; J a (l),M a ) 

M a =-J a (l) A(o) 

00 

- (a(l); J a (l), M a \ j^(-k - r«) [ dt' e^'^^' jf (k) \a(l); J a (l), M a )\E»{£' , w ) 
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where a(l); J a (l), M a is the Z-th atom ground state vector, p(X(a)) is the static weight of the 
ground state with the full set of quantum numbers A (a) except for the angular momentum, 
E a is the energy of this state, J a {l) is the total angular momentum of the atom, M a is the 
projection of the total angular momentum on the quantization axis. 

Having used the formalism developed inS for the description of decay and the lifetime of 
virtual states of quantum systems in the case, when the spectrum of the Hamiltonian of the 
system is continuous, we obtain the following formula 

oo 

(ail); Ja(l), M a \ j? (k) J dt'e-^-^-^j^i-k - r«) \a(l); J a (l), M a ) 

o 

oo 

-{a{l)- J a (l),M a \j^(-k-T^) J dl/e^'- E ' + ^j^(k)\a(iy, J a (l),M a ) = | E 

Q Jfc,Af(, 

J P (E b )dE b (o(0; Jail), M a \ if (k) \b(l); J b (l), M b ) (a(l); J a (l), M a \j^(k + r« |6(0; Ml), M b }* 

x — j- + - V / p(E b )dE b f p(E c )5(E c -hu- E a )dE c 

E a + hw-E b + R ( b +) (E a + hw) 2* jb J^ cMc J J 

x (a(l); J a [l), M a \ jf (k) \c(l); J c (l), M c ) (a(l); J a (l), M a \ j l>fl (k + r« |6(0; Ml), M b )* 
xRit\Ec) 1 — m (4.8) 

Here p(E b ) and p(E c ) are the densities of the exited atom states with energies E b and E c 
accordingly, R { b +) (E) = lim (b(l); J b (l),M b \ R(E + it) \b(l); M l ), M b) = D b (E) - U b (E) is 

the matrix diagonal element of the level shift operator determined by the equation^ 1 



R(E) = V mt + V ira ±-J£-R(E), 
E — Hi 

where Vi nt is the Hamiltonian of the interaction of the Z-th atom in the exited state with the 
crystal, P b = \b(l); Jb(/),M b ) (&(/); J&(Z),M(,| is the operator of the projection on the state 
b(l); Ml), Mb- The imaginary part of matrix element R h (E) is determined by formula 

lmR[ + \E) = ME) = y = nJ^J p(Ec) \(c(l);J e (l),M e \RM(E c ) \b(l); J h (l), M^ , 

where the summation is taken over all decay channels of the excited atom state and is the 
essence of V (the width of the crystal atom excited state). The term in the first part of 
formula (|4.8[) is proportional to R,(E a + fko) and takes into account the interaction of 



L cb 

the atom in the excited state with the crystal and in particular allows to account such an 
effect as EXAFS in the crystal polarization tensor. To describe this interaction it is enough 
to let R { c t\E c ) = (c(/); J$\Mc\ V |6(0; M l ), M b), where V is the Watson pseudo-potential 
determined by the equationu: 

v = J>|fHo>+ Y, (o\fi>^p-ff.\o) + --', 
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where Ty is the matrix of scattering by I' -th crystal atom, |0) is the crystal ground state 
vector, P = |0) (0| is the operator of the projection on the crystal ground state, - is the 

n d 

propagator (for details look at^). Here we will restrict ourself only with the first term of 
( |i.8| )-th right part since the description of the contribution of interaction of exited atom 
with the crystal to the crystal polarization tensor is worth writing a separate article. 



For the fu: 
decomposition 111 J ~ 



calculations we will use the atom current density operator multipole 



>' d ^ L ±^l -_M LM k)F^f(k) + jJ^Am^F^^)} 



(2L+1)\\^L(L + 1) w ' L ' M v ' (2L + 1)!! 

where L is the multipolarity, — L ^ M ^ L; E L>M (k), M LjM (k) and Ql,m^) are the operators 
of the electric, magnetic and charge momenta of a multipole atom; F^'jj^' (k) are the 

spherical vectors of electric, magnetic and charge typeElfl, k =— is the unit covector in the 
covector k direction. Since the field of X-ray quanta in the crystal with a high precision can 
be treated as transversal one the components with Ql,m(^) i n this decomposition can be 
neglected right away. As the magnet multipole transitions in atoms are much smaller then 
the electric multipole transitions we will content our self only with the electric multipoles. 
Having used the spherical vectors F}/ A/ (k) properties, Wigner-Eckart and Ziggert theorems, 
we obtain the equation 

— i- J2 (a(l);J a (l),M a \j?(k)\b(l);J b (l),M b ) 

a + M a ,M b 

x (o(0; J a (l), M a \j lt(t {k + r« \b(l); J b (l), M b )* 

x (a(l); J a (l)\Q L (k(\)) \b(l); J b (l))* nj(£, k(l); L) (4.9) 

where (a(l); J a (l)\ Q L (k) \b(l); J b (l)) and (a(l); J a {l)\ Q L {k{\)) \b(l); J b (l)) are the reduced 
matrix elements of the multipole charge momenta of the crystal unit cell /-th atom; 

k(l) = — —— — is the unit covector along covector (k + r^) , k = J g al3 k a k/3,k(\) = 

y g al3 (k a + Ta )(k/3 + Tp ) are these covectors' lengths; ||n^(k,k(l);L)|| is the polarization 

matrix with its elements in spiral bases e 1 = — r=(ni ± ri2),e x (l) = — =(ni(l) ± n2(l)), 

v2 v2 

where n 12 and n 12 (l) are the orthonormalized covectors in the planes orthogonal to k and 

k(l) accordingly are expressed through the Wigner functions 

L 

nJ(k,k(l);L)= Yl ^, A /(k(l))^,M(k). 

M=-L 
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If to choose the axis z in the covectors k and k + r(l) plane and to mark these covectors' 
angles relative to the axis z as $ and i?i, the formulas for the elements of the polarization 
matrix are sufficiently simplified: 



ngfr 2(1); L) = J2 ^,M(k(l))^V(k) = %Jfl - A) (4.10) 

M=-L 

By substituting the formulas (f4.8|)-( p~T0"D to the expression for j^-^w), we obtain 



IUJ i i J ^ 

tEpM0)E( (L+ / } 1 x 



2-7.(0 + - A( ,„ 

— " '■ c^fcfl)* ( a (0; J«(0I Qi(fc) [6(0; -MO) («(0; ^(01 &(*(!)) |6(0; W 



E I p(E b )dE b - 



, £ Q + ^ - # 6 + d(# + M - qr(E a + hu) 

Jb — \ L '~ Ja\ 

xd L ^(#-A)E^',cu) (4.11) 

Then the polarization tensor of the deformed crystal in the X-ray frequency range is deter- 
mined by the equation 



l i 



+ 2 



A(a) 

,3 kc , fci( e-^) c2fcLfc (i) L (a(0; Ja(0l OlW MO; jM) (a(0; J a (Q| Ql(M1)) l&(0; -MO)* 

E a + hu-E b + D(E a +Hu)- i%T(E a + ftiu) 
x4 M (tf-^)} (4.12) 



Let us discuss the obtained expression. First of all let us note that in ( |4.12|) the covector 



k and the frequency ui are not bound together by any equation. As it is seen from ( |4.7[ ), 
( [4.11] ) and ( |4.12| ), the crystal polarizability in the X-ray frequency range is a scalar value and 
the relation between the current induced in the crystal by X-ray quanta and electric vector 
or between polarization vector and electric vector is local only if we neglect the so-called 
dispersion correction and, consequently, we neglect the crystal absorption of X-ray quanta. 
If we want to account the absorption in the crystal, the statements like "the polarization 
vector in the X-ray frequency range is related to the electric vector by the equation P(r, ui) = 
x(r,u;)E(r,u;), where x(r, ui) is the crystal polarizability, and is being a scalar periodic 
function" are just incorrect. When we take into account the dispersion correction, we must 
understand that the crystal polarizability in the X-ray frequency range is in principle the 
tensor but not the scalar value and the relation between the polarization vector and the 
electric vector is unlocal. That is why the procedure of derivation of the equations ( |2.8|) and 



Q3.8| ) in parts 1 and 2 (they describe the dynamic diffraction of the electromagnetic fields on 
the media with perfect periodic and quasi-periodic structures) accounts the nonlocality of 
the response of the medium to the electromagnetic field and is also true for the X-ray quanta 
diffraction on the deformed crystals. The polarizability tensor of the medium is not being a 
member of the equations (|2.8|) and (|3.8|) but his Fourier-components Xi(k, a;), that have the 



2 ' k* 

relation between the wave covector k and the frequency u already defined: k = k^k 1 = — -. 

c 2 

If to note that electromagnetic fields with wave covectors k and k + t^' interact effectively 
only in vicinity to the Bragg conditions k 2 = k(l) 2 = (ki + t\ )(k % + r' ' l ), one can easily 
obtain the following formula for Xi(k, oj) from ( [4.121) : 



i 
, 4iTN cell 2n 



4we 2 N, 



cell 



■muj' 



fi(r {l M 



ijj- 



2J a (l) 



1 ^ 
A(o) 



p(a(o))x; 



(L+l) 



(2L + l)\\y/L{L + l 




p{E b )dE b 



J b =\L-J a \ 



2/„2L 



c 2 k 



X 



(a(l);J a (l)\Q L (u)\b(l);J b (l)) 



E a + hw-E b + D(E a + hw) - i^T(E a + hw) 



4,,^ -a)} 



It is worth to note that for 1 = the condition $ — $i = is true, and as dp JO) = 5@, 
then xi = o(k, to) is definitely a scalar value. For the wave covectors k and k + r^' (with 
magnitudes in vicinity) near the Bragg condition the •& — d\ = A9 condition takes place, 
where A6 is the angle width of the Bragg "table", i.e. is the value about 10~ -j- 10 -6 
radians. Consequently for the wave covectors k and k + t^' near the Bragg condition we 
can let dp (A9) = 5^ with high precision and can consider Xn\ „(k, u) being the scalar value 



W k ' 



w 



X(iy 



u))5p. Then for X(i){w) we have the final formula 



X(i)(?.w) = 2^ e e { 



4ne 2 N ce u 



AttN, 



cell 



1=1 

da 



mio* 



fi(r 



Oh 



(ion) 
I 



(E a + ftw-E b + D(E a + hu)) 



dE b (E a + hiu-E b + D(E a + hw)f + 



h?r 2 



■dE b + iaf on) {uj)]} 



(ion) , 



with a\ Lun> {uj) determined by the expression 



(4.13) 



a 



(ion) 



2ty 



oj 



2J a (l) 



1^ 

\(a) 



p(Ka))J2 



(L + l) 



{2L + l)\\y/L{L + l) 



\L+J a \ 



,2L-1 



(4.14) 



x J2 p(E b ) (a(l);J a (l)\Q L (u)\b(l);J b (l)) 5(E a + fku - E b + D(E a + hw))dE b 

J b =\L-J a \ J 

and being nothing but the full cross-section of the photo-ionization (photo-absorption) of 
the /-th atom in the crystal unit cell and 
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da[ on) _ 
dE b 

\L + J a \ 
J h =\L-J a \ 



^^ (A(fl)) ^l ( 2 uiSU 



P (E b ) (o(0;Ja(0IQiM|6(0;^(0) 



,2L-1 



(4.15) 



is the differential photo-ionization cross-section. These values are presently rather well 
analyzed both theoretically and experiment allyt^. It is worth to underline that while deriving 



formula ( 4.13 ) we didn't use the condition of the long-wave assumption ^a<l with a being 
the atom linear size since this condition for the X-ray radiation with the wavelength about 

o 

1 A is not true. 



V. TWO- WAVE DYNAMIC DIFFRACTION ON THE ELASTICALLY 

DEFORMED CRYSTAL. 

In this part we will discuss the X-ray quanta diffraction on the elastically deformed crys- 
tals. As the X-ray quanta field in crystals can be considered as a transverse field with high 
accuracy, the dynamic diffraction on the deformed crystals in the two-wave approximation 
according to the formulas ( |3.7|) , ( |3.8| ) and ( }4.13|) is described by the following equations: 



,('), 



E(£,w) = E ^,Lj)e ik ^ +ul ^ + E^,u;)e i[k ^ ]e+ik ^^; 

E (£,u) + vWmX(o)(£,w)Eo(£,w) + V / Rl)|X(i)(e,w)E 1 (e, 

Ei(£,w) + v¥^X(o)(£,w)Ei(£,uO 

) = 

du 



J_i (k(0 + T o»,| ) + 1 



to 



+ V / l^(OIX(-i)(e^)E (e,^ 

By definition for the elastically deformed crystals the condition 
Then owing to the fact that |x(o,i)| ~ 10 -^lO -6 



(5.1) 
<C 1 holds true. 



we can assume in ( |5TT| ) y \g\X{±\) ~ X(±i), 
\/\g(t;)\X(o) ~ X(p)- Because the deformation field and accordingly the metric tensor g^ can 
be considered to be the constant values in the area with linear dimensions about the X-ray 

1 ld(Jgg lj ) 

quanta wavelength, the constant terms containing — — - — can be neglected. The 

y/gk d£ l 

1 ■ d 2 
terms — g aH Q ^. E ,i(£, u) in the assumption of slowly varying amplitudes can also be ne- 



k 2- Q^lQ^j 

glected. If one introduces the designations n\ 



"':. 
j 



for the guiding cosines of wave covector 



k of X-ray quanta incident on the perfect crystal and n\ 



(i> 



K{ + Tj 



">i ~r Ti 



kVTT 



for the 



a 



guiding cosines of wave covector k of X-ray quanta diffracted on the perfect crystal, then 
by accounting the assumptions made above and the fact that in the Lagrange co-ordinate 

8% H — — J 



20 



as the wave covector k guiding cosines in the Lagrange co-ordinate system, the equations 
system (|5.1| ) can be sufficiently simplified 

E(£,w) = Eo(^)e* fc K ( °V+nfV(0) + ^^(^l+^^Vtt)); 



2i l n ^^9 iJ -^ + X(0)(M 



E o (e,^) + X(i)(e,^)Ei(e,^) = 0; 



{2i T [Vl + anf+n 



0) 



0) du^), « d 



,(°) 



fc 



-[VlT^n^ + 



(i) 



n 



" 9^ J ^ 

•9^ 



+ X(o)(f,w) + l 



(0) 



0- «"(£)■ 



(i) 
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+X(-i)(£,u;)E o (£,u/) = 0. 



(5.2) 



This equations system describes two-wave dynamic diffraction of X-ray quanta on the de- 
formed crystals and is the system that will be considered as basic one. 

There is Takagi -Taupin equations system that is usually used in the X-ray optics works. 
We'll show that the Takagi -Taupin equations system can be obtained from the system 



(5.2) after a set of simplifying assumptions. Let us make the following assumptions: 1) 



9ij ~ Sij + 6, 
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de 



(we recall the metric 



tensor formula gij = 8ij + 5, 



du»(Z) , . du»(0 , . du»(0du"(0 
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values of o( — ^- ^- — ) order); 
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+ s, 
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, i.e. we neglect all the 



,(°), 
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.,, 9 



0^(0, « ^ 



(0) 



dt* 



d$ 



MZM 



rif^J-E^u,); 

3) X(o,±i) (£> w ) ~ X(o,±i)(k05 i- e - we neglect the dependence of Debye- Waller factor on 
the co-ordinates. Placing these approximated equations into the system of equations ( |5.2|) , 
we obtain for the slowly varying wave amplitudes 



4^'Jj+x(°)M 



d& 



Eo(^,a;) + X(i)(w)E 1 (e,w) = 0; 



1 (i) „■,■ d 1 /- (!)„■, du u (£) 

+X(-i)(w)Eo($,w) = 0. 



a + X(o)(w) 



Ei(^w) 



Next, let us introduce vectors of 7r— and a— polarization e*?^(£), being standard for the 
X-ray optics (it is taken into account that we are applying the Lagrange co-ordinates system 



and polarization vectors are the co-ordinate functions). The values 2i— n\ S lJ T^ v (^) and 



(°), 



k 



,00 



2i— n\ \ ' S^Tj^) , where r^,(£) is the Kristoffel symbol, that appear during the polarization 

vectors differentiation and can be neglected since the metric tensor remains a constant value 
in the area with linear dimensions about the X-ray quanta wavelength. Then if we assume 
yl + a ~ 1 in the second equation, these equations look as follows 
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(i). 



+C(0x(-.)(w)^(^w) 










-a + X(o)(w) 



*i(£,w) 



(5.3) 



where C(£) is the polarization factor: Cf ff -\(£) = 1 for cr— polarization. For ir— polarization 
we have a bulky formula 
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As |x(o,i)| ~ 10 4 -T- 10 6 , we can assume C(y)(£) = n| S^rij i.e. we can assume C(y)(£) to 



be the same as that for the perfect crystal. In this case the equations (|5.3| ) for the slowly 
changing amplitudes Eq.i(£,uj) fully coincide with Takagi -Taupin equations. It is worth 
mentioning that in works mentioned above, in the formula for the X-ray quanta field in the 
crystal 



,m t 



,(°)„ 



E(£,w) = E (£,o;)e ifc(n ° (0) ^ OV(?)) + E l (^,iu)e ik{VT: ^ n ^^ +n ^" ul ^ ) 

the term n\ u l (£) in the expression for the passed and diffracted wave phases is left out 
for some reason and it is the expression that determines the X-ray quanta focusing by the 
distorted crystal both in Laue and in Bragg geometries. 

To determine the application field of simplifying assumptions made above and accordingly 
to determine the Takagi -Taupin equations application field let us determine the plane- 
parallel plate with thickness h. Let this plate to be deformed by the bending momentums 
Mi and M 2 uniformly distributed over the plate sides. Then the deformation field can be 
introduced by the following way": 
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bn(e) 2 +b 12 (e) 2 -bn(e+^) 2 +b ie ee 
+- [mo 2 + b 22 {e? - Me 3 + b 2 + b 2 ,ee 

P2 L z 



P\ 



(5.4) 
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where — = — — — ; — = — — — ; bn U = -^-; c„ u (u,v = 1, ••• ,6) is the crystal elastic 
P\ h* p 2 h s c 33 

constants. The Fig.l shows this deformed silicon plate for the parameters p\ = 2 x 10 2 cm 

and p 2 = 5 x 10 2 cm and the thickness h = 0.1cm, 
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FIGURES 




x cm 



FIG. 1. Curved Bragg mirror. 

and Fig. 2 shows the average radius of curvature of reflecting external surface of this 
plate. 
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R cm 




FIG. 2. The mirror average curvature radius 

As Fig. 2 shows, the bent silicon plate average radius of curvature with such choice of 
parameters depends on the co-ordinates of a point on the reflecting surface and is being 
the value about 3 x 10 2 cm. Let us discuss the most simple situation: the absorption of 
X-ray quanta passing through the deformed silicon plate (Bugger law). For this situation 
from equations (|5.2|) we obtain (from equations 



24^ (0) (0^(0|j + X(0)^) 



E o (£,w) = 0, 



for the simplicity we will neglect the dependence on Debye- Waller factor on the co-ordinates 
and from Takagi -Taupin equations 



*j^ 0) ^+X(o)(w) 



E o (£,u) = 0. 



Then solving the first equation by the characteristics method, i.e. by solving the system of 



the ordinary differential equations 



d? 



ds 



n 



(0) 



(0# y '(0 witn tne initial conditions ^(0) 



£o>£ 2 (0) = £o> £ 3 (0) = 0> where (£o,£o>0) are the entrance co-ordinates of the X-ray beam 
incident on the crystal, we obtain the following equation for transition function: 



T(e ,e ,h)=e 



-fcIm( X( o)(w)M£<Uo^) 



where h is the perfect crystal length. The X-ray quanta optical path s(£q,£o,/i) in the 
deformed crystal can be obtained by solving the equation £ 3 (s, £q, ^q) = h, where ^ 3 (s, ^q, ^q) 



25 



is the solution of the ordinary differential equations system mentioned above. For the second 
equation, the solution i.e. the gating function can be easily obtained and does not differ 
from the perfect crystal plane-parallel plate at all. 

rp _ e ~ fcIm (X(0)H) 51^57 

where 9 is the angle between the X-ray quanta incident on the crystal and the normal to 
the perfect crystal entrance plane directed inside the crystal. Fig. 3 show these transition 

IX 7T 

functions for 6 = — , ip = — , h = 0.01cm and quanta energies 15KeV. 
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FIG. 3. The silicon plate gating function. Here y = £ and y = £, ■ 

At the upper surface of this figure is presented the transition function for the plane- 
parallel perfect silicon crystal plate and the lower surface describes the X-ray quanta ab- 
sorption by the deformed crystal. The distinction shows itself in the term of the second order 
and is very sufficient for such a thin plate. It means that even in such a simple situation as 
the absorption of X-ray quanta passing through the deformed crystal, the Takagi -Taupin 
equations don't correctly describe the real experimental situation (this distinction can be 
neglected only for the crystal deformation with the average radius of curvature more then 
40m). 

To finally define the application range of Takagi -Taupin equations let us compare two 
terms: 

1) 1 - [VT+amP + nf^^-]g lj [VlT^nf + n^S in the second equation of 
the equation system (|5.2| ) and 

2) —n\ S^Tfj, a in the second equation of the Takagi -Taupin equations system 

( |5.3j). Figure 4 shows the difference 
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'" -l ,^ u "' ^U M^Tzrr^f + n y) ^l ] _ L n f#. (// 



- 1 - I n (D^ T a*"® 



between these values for the deformed silicon plate shown on the Fig.4-6 for the energies 
5KeV, lOKeV and 15KeV of the quanta. 



5 KeV 



10 4 D 




y cm 



FIG. 4. The value D for the quanta with energy 5KeV . Here y = £ and y = £ . 



10 KeV 




y cm 



FIG. 5. The value D for the quanta with energy lOKeV. Here y = £ and y = £ • 
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15 KeV 



10 5 D 




FIG. 6. The value D for the quanta with energy 15KeV. Here y = £ and y = £ . 

If one notify that for energies 5KeV,10KeV and 15KeV we have Re(x(o)( w )) = 
-0.394922 x 10~ 4 , -0.976232 x 10~ 5 and -0.431670 x 10~ 5 accordingly, this figure shows 
that by accepting the simplifying assumptions in the equations (|5.2|) leading to the Takagi 



-Taupin equations we reject the values that are 10 times greater than the values we keep 
(for the crystals with the average radius of curvature less than 10m) while obtaining these 
equations. This difference shown in Fig. 4 becomes less then Re(x(o)( a; )) f° r the deformations 
with the curvature average radius more then 40m, i.e. very slightly deformed crystals. That 
is why one should use exactly the equations system ( |5.2| ) for the theoretical analysis. This 
system is a little bit more complicated but is suitable for the description of the X-ray optics 
devices being met with in practice. 

VI. BEND BRAGG MIRROR. 

The equations Q5.2p ( equations for the slowly varying amplitudes of the field of X-ray 
quanta in the two-wave dynamic diffraction on the deformed crystals) for the arbitrary 
deformation field can be solved in general only numerically. In the article we underlined 
more than once that the X-ray optics can be considered as the geometrical optics from 
the point of view of general electrodynamics. That is why to obtain the coefficient of Bragg 
reflection from the distorted Bragg mirror we will use the geometrical optics locality principle 
(see for ex.E-3) that evidently was first time introduced in electrodynamics by V.A. Fockp3 
who investigated the reflection of the electromagnetic wave with smooth arbitrary shape 
front on the smooth arbitrary shaped surface. According to the locality principle the wave 
reflection in the every mirror point acts like if the incident wave is being smooth one and 
the curvilinear surface in this point is replaced by the tangent plane. But it is necessary to 
take into account that the incident and the reflection angle change their values from point 
to point on this surface. Consequently the locality principle allows us to use the reflection 
coefficients obtained for the plane waves and plane surfaces, and in particular the Bragg 
reflection coefficient for the diffraction X-ray optics as the primary approximation. The 
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Bragg reflection coefficient in the case of the plane wave reflecting from the parallel-sided 
perfect crystal thick plate is determined by formula: 



n 



Cx(-i)M 



a 



x (0 )M(i + /?)--± 



X(o) 



uj 



(6.1) 



0) 



CWM 



where C is the polarization factor mentioned above, (3 



n 



(i) 



n 



(0) 



is the asymmetry parameter, 



-k 2 



k 2 



is the parameter of precise Bragg condition 



X 2 M = X(-i)(w)X(i)M, a = - 
deviation. 

Let the X-ray quanta wave packet incident on the Bragg mirror with arbitrary smooth 
deformation field u(£) has the electric vector for quanta with frequency uj ( or the energy e 
) determined in the Lagrange co-ordinates system by the equation 

where n\ are wave vector guiding cosines and E'°^(£, a;) and ^fo(^,s) are the slowly chang- 
ing amplitude and phase of the X-ray quanta incident on the mirror. Then according to 
the locality principle of geometrical optics this reflection coefficient of the wave packet for 
distorted Bragg mirror has the functional form like that for the plane wave and plane mirror. 
But it is necessary to take into account that the values X(o,±i) , C, [3 and a for the distorted 
Bragg mirror change from point to point on this mirror surface. If the dependence X(o,±i) 
on co-ordinates is determined by formula (|4.13|) , the dependencies of C,/3 and a on the 
co-ordinates can be defined by the following ideas. The wave covector components of the 
X-quanta incident on the distorted Bragg mirror in the Lagrange co-ordinate system to be 
determined by the equation 



m\e) = k 






nf ] +nf 



+ 



e=o 






C 3 =o 



Then the wave covector X-ray quanta components in the Lagrange co-ordinates system 
having diffracted on this mirror are determined by formula 



m\n + n = k 



VW + *.(5fi) + 



d*o(M 



d? / C 3 =0J 
According to the definition we have for parameters ck(£ ,£ 2 ),/3(£ ,£ 2 ) and C(^ ,^ 2 ) 
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where n\ '(£ , £ ) and wh (£ , £ ) are the wave covectors guiding cosines of X-ray quanta 
falling and diffracted on the distorted Bragg mirror surface in the Lagrange co-ordinates 
system: 
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Consequently the Bragg reflection coefficient of the distorted mirror in this approximation 
looks like 



1 £2\ 



ne,z 



c X (-i)(e,e;v) 



X( o)(e, e; w)(i + p) - ^p- ± A(e, a 



1 t 2\ 



A(£\£ 



1 c2. 



X(o)(^,^ 



^(1 + /?)-^ 



ewe* 1 ,? 2 ;") 



(6.3) 



Here X(o,±i)(£\£ 2 ;w) = (x (0 ,±i)(f; w)) $3=0 , X 2 (£\£V) = X®(t\?',u)X{.-V)(?,?'M- The 
asymmetry parameter /3 and the polarization factor defined by formulas (|6.2| ) can be con- 
sidered the same as for the plane mirror in the first approximation. If besides we neglect the 
Debye- Waller dependency on the co-ordinates, the formula ( p.3| ) for the Bragg coefficient is 
extremely simplified: 
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CX ^ U) (6.4) 



The fact, that the distorted Bragg mirror with rather big average radius of curvature can be 
represented as a set of plane Bragg mirrors turned relatively to each other, is widely used 
in X-ray optics systems designing. Nevertheless we haven't found the formulas (|6.3[ - |6.4|) 
introducing this representation (or their analogues) in other papers. The linear dimensions 
of the area participating in the forming of the diffracted field in the point (£\£ 2 ) on the 
crystal surface are about the X-ray quanta crystal absorption length l ahs ~ 10~ 2 -j- lCT 3 cm. 
If the deformation field in this area can be considered as the constant value, the formulas 



( |Q| - p.4D are quite applicable. 

Once the Bragg reflection coefficient is obtained, we can define the field of X-ray quanta 
diffracted on the distorted Bragg mirror reflecting surface 



The equations ( |6.3| - |6.4|) and ( |6.5|) solve the problem of X-ray quanta two-wave dynamic 
diffraction on the elastically deformed crystal in the Bragg geometry, i.e. on the incurved 
Bragg mirror. The only problem remaining is the X-ray quanta propagation and focusing 
in vacuum (or in the air), i.e. the Hemholts equation solution. 

(A + k 2 )E(r,uj) = 

with the assumption that on the distorted Bragg mirror surface there is the X-ray quanta 
field A {0 \r,u)e lkSo{r) defined with the amplitude A {0) (r,uj) = TZ^\^)E {0) ^\^;uj) and 
phase S (r,e) = $(i)(^ 1 ,^ 2 ,£:). Beyond the caustics this task can be solved by geometrical 
optics methodsE-j leading to the formula 



E{r,u>) = Hi?, ?)&(?,?;«;) ,/ (g _ ^ _ — e^+ikS^e) (6 . 6) 



where 



s(r,t\t 2 ) = VW 1 - 1 1 - uHt\e, o)) 2 + (* 2 - e - ^W, ¥, ow + p - u*(t\?, ow 

(6.7) 



is the distance from the Bragg mirror surface point with parameters (£ , ^ 2 ) to the point with 
co-ordinates (a; 1 , x 2 , x 3 ), and R\ and i?2 are the curvature radiuses of the wave diffracted on 
the Bragg mirror surface defined by the equations 
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(6.8) 



Here # u is the fully asymmetrical tensor of the third rank (<5 123 = 1), iVj(£ , £ 2 ) is the 
normal vector to the X-ray quanta wave front on the Bragg mirror surface in the Euler 
co-ordinates system defined by the formula 
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where J (£ , £ ) is the matrix inverse to the Jacob matrix J(£ , £ 
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(6.9) 



§ 3 =0 



and Xj(£) = £ + w'(0 are the Euler co-ordinates of the deformed crystal points. 

In the points where the equations S = R12 are fulfilled, the geometric-optical solution 
(40) turns into the infinity (has a singularity). These points are called the focal ones and the 
set of the focal points is called the caustic. Consequently the multitude called the caustic is 
determined by the equations 



.1" 



(/)1.2 



(e + « j (o) f 3 =0 + r^z 1 , ?)***(?, e 



(6.10) 



As the example illustrating formulas ( O )-( BTT^ ) we will consider a simple situation when 
the plane waves packet of X-ray quanta with parallel vectors and energies in the interval 
from HKeV to YJKeV incident on the plane-parallel silicon plate with thickness h = 0.1cm 
and reflection vector r(l, 1, 1) perpendicular to the plate entrance surface. We will consider 
the crystal plate oriented in such the way that the Bragg condition precisely holds for the 
X-ray quanta with energies YhKeV . Next this plate is distorted by momentums Mi and 
M 2 uniformly distributed over the plate sides in such a way that the deformations field is 



determined by the formulas ( |5.4j ) from the previous part and the distorted Bragg mirror is 
the result of this deformation shown on Fig.l. Since by accepted condition the packet of 
plane waves is incident on this mirror, one would let \&o(£, £) = in formulas ( |6.3| )-( |6"TT0"D 
in this situation. The angle between the diffraction plate and the co-ordinate plate (£ 1 ,£ 3 ) 
will be denoted as (p. Fig. 7-10 shows the graphs obtained with the help of formula ( |6.4|) 
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for the module of coefficient of Bragg reflection from this distorted mirror for X-ray quanta 
energies llKeV, 13KeV, 15KeV and 17 KeV with <p = 0. 
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FIG. 7. Bragg reflection coefficient. Here y = £ and y = £ . 
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FIG. 8. Bragg reflection coefficient. Here y = £ and y = £ . 
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FIG. 9. Bragg reflection coefficient. Here y = £ and y = £ . 
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3.0 62 
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FIG. 10. Bragg reflection coefficient. Here y = £ and y = £ 



We remind that the caustics (the focal point multitude) of X-ray quanta diffracted on the 
Bragg mirror is defined exclusively by the X-ray radiation field and does not depend on the 
reflected quanta intensity distribution (Bragg radiation coefficient) on the crystal surface. 
By Fig. 11-14 we give the caustics for X-ray quanta energies HKeV, 13KeV, 15KeV and 
17 KeV obtained from formulas ( |6.9[ - |6.10| ) for the whole Bragg mirror (without accounting 
the intensity distribution of quanta reflected on the whole surface of the mirror, i.e. without 
accounting the Bragg reflection coefficient dependency on the co-ordinates of the point on 
the mirror surface). 
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FIG. 11. Diffracted quanta caustic. 
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FIG. 12. Diffracted quanta caustic. 
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FIG. 13. Diffracted quanta caustic. 
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FIG. 14. Diffracted quanta caustic. 



It is clear that in experiment there will be presented only the part of this caustic surface 
which X-ray quanta are coming at, once reflected from the mirror area with the Bragg 
reflection coefficient sufficiently differing from null. These areas for our example are shown 
on fig. 5. Fig. 15-18 shows caustics (focal points manifold) obtained also with the help of 
formulas ( |6.8[ - |6.10|) but taking into account the dependency of Bragg reflection coefficients 
on the co-ordinates. 
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FIG. 15. Caustic at 9.46 < f 1 < 9.475. 
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FIG. 16. Caustic at 3.994 < f 1 < 4.013. 
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FIG. 17. Caustic at -0.007 < f 1 < 0.004. 
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FIG. 18. Caustic at -3.069 < ^ < -3.06. 

It is worth mentioning that even in the simplest case we deal with rather compli- 
cated (that is clear from these diagrams) combination of elementary differentiable peculiar- 
ities (catastrophes) Ai, A2, A$, A± and D 4 (by V.I. Arnold classification) of phase function 
$(i)(^ 1 , £ 2 , e) + S(v, £ , £ 2 ). That is why one can not even talk about the fit (frequently met 
in the X-ray optics articles, see for ex.") of X-ray quanta intensity distribution near the 
caustics (focus) with Gaussian functions . The caustics intensity distribution is described 
by special catastrophes functions the simplest of them being Airy functions and Piersey 
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integral" investigated in details in articles^. If one models the intensity distribution by 
Gaussian beams, it is necessary to use the special method of Gaussian beams summing 
discussed in". 

VII. ACKNOWLEDGMENTS 

We are indebted to Prof. S. Maksimenko and Prof. G. Slepyan for stimulating discus- 
sions. 



39 



REFERENCES 

x Takagi S., Acta Cryst. 15, 1311-1312, 1962. 

2 Tuipin D., Bull. Soc. Fr. Mineral. Cristtalogr. 87, 469-511, 1969. 

3 Chukhovskiy F. N., Gabrielyan K. T., Petrashen P. V, Acta Cryst. A34, 610-621, 1978. 

4 Chukhovskiy F. N.,Malgrange C, Gronkowski J, Acta Cryst. A52, 47-55, 1996. 

5 N. Kato, Acta Cryst. A46, 672-681, 1990; Acta Cryst. A48, 829 -841, 1992; Acta Cryst. 

A 203-213, 1997. 
6 A.S. Ilyinsky, L.V. Keldysh, Electromagnetic field interaction with the matter, Moscow 

Univ. Publ.,M.,1987 
7 N.A. Hizhnyak, Integral equations of macroscopic electrodynamics , Navukova Dumka, 

Kiev, 1986. 
8 Bjorken J.D., Drell S.D., Relativistic quantum mechanics, McGraw-Hill Company. 
9 Goldberger M. L., Watson K.M. , Collision Theory, John Willey&Sons, Inc. 1964. 

10 G.L. Birr, G.E. Picus, Symmetry and deformation effects in semiconductors, Nauka, 
Moscow, 1972. 

11 Bohr A., Mottelson B. R., Nuclear Structure, vol. 1, W.A. Benjamin Inc., New York, 
Amsterdam, 1969. 

12 V.M. Dubovik, A. A. Cheshkov, Physics of Elementary Particles and Atomic Nuclei, 1974, 
vol.5, p.791. 

13 D. A. Varshalovitch, A.N. Moscalyov, V.K. Hersonskiy, The quantum theory of angle mo- 
mentum , Nauka, 1975. 

14 Rose M.E. Multipole fields, McGraw-Hill Company, London, 1955. 

15 Amusia M.Ya. Atom photoeffect, Nauka, Moscow, 1987; Amusia M. Ya, Kheifets A. S. , 
Phys. Lett., 1981, v.82A, p. 437. 

16 S.G. Lekhnitskiy, The elasticity theory of anisotropic body , Nauka, M.,1979. 

17 Yu.A. Kraftzov, Yu.I. Orlov, The geometry optics of non-uniform media, Nauka, 1980. 

18 V.A. Fock, Electromagnetic fields diffraction and propagarion problems, Sovietskoie radio, 
1970. 

19 V.I. Arnold, A.N. Varchenko, S.M. Gussein-Zade, The peculiarities of differetiable maps, 
Nauka, M., 1982. 

20 V.I. Arnold, A.N. Varchenko, S.M. Gussein-Zade, The peculiarities of differetiable maps, 
Nauka, M., 1984. 

21 Chukhovskiy F.N., Chang W. Z., Foster E. , J. Appl. Phys. 77(4), p. 1843-1848. 
22 Pearcey T., Philos. Mag. 1946, 37, p. 311-317. 

23 A. S. Kryukovskiy, D.S. Lukin, E.A. Palkin, Uniform asymptotics of quickly oscillating 
functions with singular saddle points, A.Sc. USSR, Radiotechnics and Electronics Insti- 
tute, Preprint # 41 (413), M., 1984 

24 Grikurov V. E., Popov M. M., Wave motion, 1983, v. 5, p. 1181-1191. 



40 



